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We investigate the nature of the dynamo bifurcation in a configuration applicable to 
the Earth's liquid outer core, i.e. in a rotating spherical shell with thermally driven 
motions. We show that the nature of the bifurcation, which can be either supercritical 
or subcritical or even take the form of isola (or detached lobes) strongly depends on the 
parameters. This dependence is described in a range of parameters numerically accessible 
(which unfortunately remains remote from geophysical application), and we show how 
the magnetic Prandtl number and the Ekman number control these transitions. 
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1. Introduction 

The origin of the Earth magnetic field remains a challenging issue for physicists. 
The accepted theory, known as dynamo theory, describes the transfer from kinetic 
to magnetic energy as an instability process. Above a given threshold electrical 
currents, and thus magnetic fields, are amplified by a turbulent flow of a conducting 
fluid. 

This problem can be described in its simpler form by a set of coupled partial 
difii'erential equations (see section 2). The state of this system is fully characterised 
by four independent controlling parameters. The Ekman number E, which can be 
interpreted as measuring the ratio of the length of the day to a typical viscous 
timescale, this number is extremely small in the case of the Earth's core (the Earth 
is in rapid rotation). The magnetic Prandtl number Pm, measuring the ratio of 
a typical ohmic timescale to viscous timescale, it is a characteristic of the fluid 
and is minute for all liquid metals (including liquid iron in the Earth's core). The 
hydrodynamic Prandtl number Pr, also characterizing the fluid and which compares 
a thermal timescale to a viscous timescale, this number is of order one, or slightly 
less, but non- vanishing. Finally the Rayleigh number Ra, which measures a ratio 
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of driving forced to forced slowing down the motion. Its value is difficult to assess 
in a simple Boussinesq model. This will be the most obvious controlling parameter, 
which needs to be varied to investigate dynamo properties. 

Kutzner and Christensen [1] performed a first systematic exploration of the pa- 
rameter space available to numerical models of the geodynamo. They produced a 
"phase diagram" for dynamo action in terms of Ekman number E and magnetic 
Prandtl number (or equivalently here Roberts number, q — Pm/Pr) and Rayleigh 
number Ra. Their study was later extended to lower values of E and Pm in [2\ . We 
present a schematic sketch of their results in figure [1] (See also the PhD thesis of 
Carsten Kutzner 3 ). Each diagram corresponds to a given Ekman number, rang- 
ing here from 2.10"^ to 10""*. Kutzner and Christensen identified three "phases" 
on these diagrams. In part of the parameter space, no dynamo solutions were ob- 
tained. In the complementary region, in which dynamo is obtained, two regimes are 
identified. One at low forcing, characterised by strongly dipolar, but non-reversing 
dynamos, and another corresponding to reversing, but multipolar dynamos. This, of 
course, is not a very good news when comparing these models with the geodynamo 
(characterized by a reversing strongly dipolar field). 

As pointed above, all models presented in these diagrams suffer from parameters 
extremely remote from their geophysically relevant estimates. In particular, the 
Ekman number E is over-estimated by a factor in excess of a billion and the magnetic 
Prandtl number Pm by a factor close to a million. In fact this last number offers a 
ratio of diffusivity which is equal to unity or larger in the numerics, whereas it should 
be as small as 10^^ in the Earth. The diagram reveals that if the magnetic Prandtl 
number is decreased in numerical models (all other parameters being fixed), dynamo 
action is lost! This does not come as a surprise, as it corresponds to increasing the 
ohmic resistivity, but appears to prevent any direct comparison of the models with 
the actual geodynamo. 

In fact the main result derived from these phase diagrams is much more op- 
timistic. The authors noted that the threshold value Pm^, below which dipolar 
dynamos are lost decreases with the Ekman number. They proposed [4l[2] a scaling 
in the form Pm^ cx E" and proposed a numerical fit to a ~ 3/4 (see [5] for a the- 
oretical explanation of this behaviour leading to a ~ 2/3). Which could indicate a 
connection in the parameter space between the models presently available and the 
regime relevant to the geodynamo. 

Part of the Kutzner & Christensen results have also been discussed in [4j. In 
this earlier study, the authors proposed that for E = 10"'^ the B = state is always 
stable. The corresponding bifurcations are therefore interpreted as subcritical by the 
authors. For E = lO^'*, the authors identify only supercritical bifurcations (because 
they find B = to be unstable). They however did not perform a systematic study 
to characterise the dynamo bifurcation, nor determined how it is controlled by the 
governing parameters. This is the object of the work reported here. 

We present here a review of the results obtained and first described in french 
in Morin 2005 [6^ (available online). These results were subsequently presented at 
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Fig. 1. "Phase diagram" for dynamo action for decreasing values of the Ekman number (after 
Kutzner, 2003). Three domains are identified. Domain I corresponds to parameters for which no 
dynamo solutions were obtained. Domain II corresponds to dynamo solutions dominated by a 
strong axial dipole. Domain III corresponds to reversing multipolar dynamos. 



several international conferences [7]. They were not published in english so far, but 
have probably inspired more recent articles [S]. We should finally stress that while 
our results were originally obtained using only the numerical code call "Magic" 
(developed by G. Glatzmaier, U. Christensen and J. Wicht), we have since then re- 
produced these results using the local "PaRoDy"code. Simulations with this second 
code have confirmed the results presented here. 
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2. Numerical investigations 

In all the simulations presented here, we solve the MHD equations in the Boussinesq 
approximation in a rotating frame between two concentric spheres (ro/r^ — 0.35). 
The reference frame is set such that the velocity vanishes on both spheres (no-slip 
boundaries), a temperature difference is maintained across the shell, and both the 
inner and the outer domain are assumed to be electrically insulating. The equations 
governing the solenoidal velocity u and magnetic B fields and the temperature T 
can be written in non-dimensional form 

E [dtu + (u • V)u] = - Vtt + E Au - 26^ x u 

+ RaTr + Pm"i (V X B) X B, (1) 

9tB = V X (u X B) + Pm^i AB , (2) 

dtT+{u-V)T = Pr-^AT, (3) 

using the previously introduced Ekman number E = ly/flL'^ , the (modified) 
Rayleigh number Ra = agATL/M , the Prandtl number Pr = u/k , and the mag- 
netic Prandtl number Pm = v/rj , with L = Tq — r.i, and g ~ go/ro where go is the 
gravity at r = Tq. For simplicity, we use here the same non-dimensional form as 
in [9]. The rapid rotation of the reference frame is measured by the smallness of E. 
In the following, the Prandtl number is set to unity, and thus q = n/rj ^ Pm. 

As we will present our results using Ra/Rac as controlling parameter, where 
Rac is the critical Rayleigh number for the onset of thermal convection, we want to 
recall, for completeness, that the values of Rac are respectively 55.9, 60.8, and 69.7 
for the Ekman numbers considered here i.e. 10~'^, 3 • 10^^, and 10^^ . 

2.1. Influence of the magnetic Prandtl number on the dynamo 
bifurcation 

Let us first fix the Ekman number and investigate the nature of the dynamo bifurca- 
tion varying the magnetic Prandtl number. With an Ekman number of E = 3.10"'*, 
when the magnetic Prandtl number is set to Pm = 6, the bifurcation diagram 
obtained numerically is reported on figure [21 We use filled dots to indicate stable 
solutions, and opened dots to indicate unstable solutions. The diagram reports mag- 
netic energies, for varying values of the normalised Rayleigh number, dots indicate 
the time averaged values of the magnetic energy, while the standard deviation is 
indicated using error bars. The simulations are initialised with a small temperature 
perturbation, a vanishing velocity, and an infinitesimal magnetic perturbation. 

One can note on figure [2] that the first point with non-zero magnetic field is 
obtained for Ra = 2.57 x Rac, the previous point on the diagram Ra = 2.42 x Rac, 
corresponds to a stable B = state. Figure[2la demonstrates how the time averaged 
magnetic energy continuously tends to zero as he controlling parameter is decreased. 
This demonstrates the supercritical nature of this bifurcation. 
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Fig. 2. Mean magnetic energy as a function of the Rayleigh number. The bifurcation is found to 
be supercritical for E = 3.10~* and Pm = 6. Error bars indicate the standard deviation of energy 
fluctuations around the mean value. Graph b is an enlarged version of graph a, centered on the 
dynamo bifurcation. 



The onset of dynamo action is here obtained close to the onset of convection. At 
the onset, the dynamo reaches a regime with constant kinetic and magnetic energy. 

We should stress that each point on such diagram, corresponds to simulations 
of a few 100 of hours on a parallel computer. Computations are three-dimensional 
and need to be carried over a long enough time to obtain well converged averaged 
values. 
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Fig. 3. Mean magnetic energy as a function of the Rayleigh number. This subcritical bifurcation 
was obtained for E = 3.10~* and Pm = 3. Error bars indicate the standard deviation of the 
magnetic energy fluctuations around its mean value (dots). A close-up on the dynamo threshold 
is presented in graph b. 



When the magnetic Prandtl number is decreased to Pm = 3, the onset of dy- 
namo action occurs between Ra = 3.85 x Rac and Ra — i x Rac (further remote 
from the onset of convection). The magnetic energy reported on figure [3] exhibits 
an hysteretic behaviour characteristic of a subcritical bifurcation. The lower point 
on the subcritical branch is "metastable" . It corresponds to Ra = 3.16 x Rac and is 
designated on figure [3] by a stripped dot. By metastable, we mean a self-sustained 
dynamo solution, which could be maintained for a long time compared to viscous 
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and ohmic diffusion times, but which eventually collapses when the system expe- 
riences a sufficiently large fluctuation (in the sense of a deviation from the mean 
magnetic energy). Once the dynamo solution is lost, the noise being multiplicative, 
the system remains in this state. In the case of this example, the dynamo solution 
was maintained for 40 viscous times, i.e. approximately 17 magnetic decay times 
before coUapsing. 




Fig. 4. Growthratc of the magnetic energy of a perturbation versus normalised Rayleigfi number 
for E = 3.10~* and Pm = 3. The growthrate of the solution first increases as the value of 
the Rayleigh number is increased above the dynamo threshold. It reaches a maximum value, 
and then decreases to eventually become negative. The hydrodynamic solution restabilises above 
Ra ~ 6 X Rac . 



An intriguing characteristic of this bifurcation diagram is that the state B = 
restabilises for large enough values of the Rayleigh number, above the dynamo 
threshold. This is not due to an isolated point which behaves unexpectedly. In 
fact one can compute the linear growthrate of the magnetic field before reaching 
saturation, and the decay rate when B = is stable (see figure |4|). This growthrate 
first increases above the onset, but then reaches a maximum as the Rayleigh number 
increases and eventually becomes negative again. 

A third type of bifurcation is obtained for a magnetic Prandtl number of Pm = 
1.5 (see figure [Ha). This bifurcation has an atypical shape sometimes referred to 
as "isola" . In this regime, the B = solution does not destabilise, but a detached 
branch of dynamo solution exists for Ra/Rac G [4.44,9.18] 

We have observed in this section, decreasing the magnetic Prandtl number from 
Pm = 6 to Pm = 1.5, that we could describe supercritical, subcritical and isola 
diagrams for dynamo action. We shall note that the Rayleigh number had to be 
increased in this process (as the effects of Ohmic diffusion increase for decreasing 
values of the magnetic Prandtl number). 

2.2. Influence of the Ekman number on the dynamo bifurcation 

Rather than studying the parameter space at fixed Ekman number and varying the 
magnetic Prandtl number, we will now vary the Ekman number for a fixed magnetic 
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Fig. 5. Mean magnetic energy versus Rayleigh numbers. Isola diagram were obtained for E = 
3.10~^ and Pm = 1.5 (left plot) and E = 10~^ and Pm = 3 (right plot). Black dots (resp. stripped 
dots) indicate stable (resp. metastable) states. Error bars indicate the standard deviation of the 
magnetic energy fluctuations around its mean value. Metastable states are unstable to sufficiently 
large fluctuations of the velocity field. The probability for these fluctuations to occur are non zero. 
Purely hydrodynamic solutions were found to be always stable for this regime of parameters. 



Prandtl number Pm — 3. Results obtained for E — 10~^ are reported on figure[5]b. 
This diagram corresponds again to an isola. All dynamo solutions were found to 
be metastable. The extreme points on this branch were maintained for roughly 6 
other points were maintained in time for up to 140 Tj^. We have seen previously 
that the dynamo bifurcation is subcritical for an Ekman number E = 3.10"^ and 
a magnetic Prandtl number Pm ~ 3. Let us now turn to E = 10^^, the resulting 
bifurcation diagram is presented on figure |6la and corresponds to a supercritical 
bifurcation. 
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Fig. 6. Mean magnetic energy versus Rayleigh number. These supercritical (a) and subcritical (b) 
bifurcations were obtained respectively for E = 10~* and Pm = 3 and E = 10~* and Pm = 0.67. 
Error bars indicate the standard deviation of the magnetic energy fluctuations around its mean 
value. 



By increasing the Ekman number at fixed magnetic Prandtl number, we there- 
fore successively described supercritical, subcritical and isola bifurcation diagrams. 
We further extended our study to E = 10""* and Pm = 0.67 (figure [6lb) and ob- 
tained a subcritical bifurcation. 



June 8, 2010 0:46 WSPC/INSTRUCTION FILE MorinDormy 



8 Morin & Dormy 

2.3. Interpretation 

We have studied different bifurcations obtained for Eknian number values ranging 
from 10^^ to 10^^, magnetic Prandtl number values from 0.67 to 6 and Rayleigh 
number values from Ra ~ 2 x Rac to Ra ~ 10 x Rac. In this parameter regime, the 
evolution of the bifurcation diagram is summarised in figure [T] For a given Ekman 
number, a supercritical bifurcation (top plot) is obtained for a sufficiently high value 
of the magnetic Prandtl number. By decreasing its value, it is possible to obtain a 
subcritical bifurcation (middle plot), which may exhibit unusual features, such as 
re-stabilization of the purely hydrodynamical state. An unstable branch therefore 
must exist, it could be connected, for larger values of the Rayleigh number, to the 
stable branch corresponding to dynamo solutions. If the magnetic Prandtl number is 
further decreased, the range of Rayleigh number for which the non dynamo solution 
is unstable vanishes. An isola is then obtained (bottom plot), in this situation the 
purely hydrodynamical solution is always stable. The very same sequence is obtained 
by increasing E at fixed Pm, as the dipolar domain (domain II on figure 1) shifts 
towards lower values of Pm as E is decreased. 

2.4. Feedback and coupling 

Let us now turn to the effects of the saturated magnetic field on the fiow. We 
compare on figure [S] simulations obtained for E = 1.10~^ and Pm = 3, before and 
after saturation. The magnetic field clearly acts here to inhibit thermal convection. 
Not only does it lower the averaged value of the kinetic energy, but it also reduces 
the amplitude of fiuctuations. This behavior is typical of the parameter space we 
investigated, there are however a few noticeable exceptions. 

In the case of the subcritical bifurcation obtained for E = 3. 10^"* and Pm = 3, 
we can investigate a range of Rayleigh numbers from Ra — 3.16 to 3.30 x Rac 
(which are lower than the linear threshold value). It is then found, figure |9l that the 
presence of the magnetic field yields stronger velocity fiuctuations. Indeed the purely 
hydrodynamical solution (without magnetic field) is presented on the right of each 
graph (either after a loss of dynamo action in the metastable case Ra = 3.16 x Rac, 
or through an independent non-magnetic simulation). 

Dynamo solutions corresponding to the isola branch described for E — 10^^ and 
Pm = 3 are metastable. Magnetic energy fiuctuations are here comparable to the 
mean value, which leads to a loss of dynamo action after a time which depends on the 
Rayleigh number as well as on the particular choice of initial conditions. We should 
highlight a simulation performed with a Rayleigh number Ra = 3 . 76 x Rac , presented 
on figure [TOl Dynamo action was here maintained for 140 (with — L"^ /v) before 
being lost, this corresponds to 20 (with = rl/rj), and to ~ 200 dipole decay 
time {Td = Tjj/tt^). Such behaviour raises questions as to "how long should a dynamo 
be investigated to assess its stability?". Once lost, the dynamo solution cannot be 
recovered, as the purely hydrodynamic solution is stable for this parameter regime. 
Once in the basin of attraction of the B = state, the field cannot recover. This is 



June 8, 2010 0:46 WSPC/INSTRUCTION FILE MorinDormy 



Dynamo bifurcation 9 




Magnetic Energy 




i i i i 



Ra/Rac 

Fig. 7. Summary of the typical bifurcations obtained numerically. As the magnetic Prandtl num- 
ber is decreased at fixed Ekman number one obtains successively supercritical (top), subcritical 
(middle) with restabilisation, and isola (bottom). The same sequence is achieved by increasing the 
Ekman number at fixed Rossby number. The dashed lines indicating unstable branches are here 
speculative, except for the B = state. 



due to the multiplicative character of hydrodynamic fluctuations in the induction 
equation. The effect of the hydrodynamic fluctuations on the field vanishes near 
B = 0. 

Let us now turn to the statistical properties of the fluctuations of kinetic and 
magnetic energies. Figure [12] reports the probability density functions (PDF) of 
both energies for Pm = 3 and on the top row E = 3.10"^, Ra = 3.7 x Rac (sub- 
critical bifurcation), on the bottom row E — 10^*, Ra — 3.3 x Rac (supercritical 
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Fig. 8. Magnetic energy (solid lines) and kinetic (dashed lines) as a function of time for E 
1.10"", Pm = 3 and from left to right and top to bottom Ra = 2.87, 3.01, 3.30 and 4.30 X Rac. 



bifurcation). In both cases the PDFs of the magnetic energy are much less symmet- 
rical than their counterpart on the kinetic energy. This is a robust characteristic, 
which was observed in all our simulations. 

This appears again to be a direct consequence of the multiplicative character 
of the hydrodynamic fluctuations in the induction equation. Whereas temperature 
fluctuations act as additive noise on the velocity equation, velocity fluctuations act 
as a multiplicative noise on the magnetic field (i.e. a term proportional to B). One 
can show on a very simple system how a multiplicative white noise (with symmetric 
PDF) yields a chi-2 law (asymmetric PDF) for the variable affected by the noise. 
Let us consider the temporal evolution of a scalar, for a conservative system, we 
can write 

dX dG 

l^ = -dx- 

An additive noise would drive fluctuations of X within a fixed potential G, a mul- 
tiplicative noise will instead modify the shape of the potential. Let us consider, for 
example, the case of a steady subcritical bifurcation in the presence of a multiplica- 
tive noise. The normal form (assuming the X ^ —X invariance) yields a potential 
of the form 



G{X) 



X' 



a 



X' 



1 



■X' 



(5) 



2 4 6 

In the absence of noise (^ = 0), solutions to (jUlS]) are: the trivial X = (stable for 
H < unstable otherwise), and the four roots oi ^ + aX"^ — X"', for —a^ /A ^ /i < 0. 
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Fig. 9. Magnetic (solid) and kinetic (dashed) energies as a function of time for E = 3.10"* and 
Pm = 3 and from top to bottom and left to right Ra = 3.16,3.30,3.70 et 3.85 X Rac For all 
graphs, the signal on the far right part of the graph corresponds to the kinetic energy evolution 
in the absence of magnetic field (purely hydrodynamic solution). 
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Fig. 10. Magnetic (and kinetic on the right graph in dashed lines) energies versus time for E = 
10"'^, Pm = 3 and Ra = 210. The dynamo solution is lost after some 140 r^. Both the mean and 
the fluctuating kinetic energies increase after the magnetic field collapse. 

Two of these are unstable, the two others being stable (subcritical branch) and 
continuous for ^ ^ . 

Let us now consider the subcritical branch of this system with noise. Below 
the threshold, the parameter /i is negative. Figure [TT] examplifies the shape of the 
potential G{X) for different values of £,{t), and thus of the prefactor of the term. 
When it is low enough (solid line), there exists a negative minimum. This minimum 
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provides the preferred state of the system. When £^{t) becomes large enough, the 
potential can fluctuate until both minima of the potential reach the same value 
(dotted curve). If the prefactor of the term becomes even larger, the potential 
minimum associated to non-zero X becomes local and the prefered state of the 
system becomes X — 0. 

It is important to appreciate that X = is an absorbing state. Indeed the noise 
being multiplicative will not affect the system any longer when X = O.li fluctuations 
are large enough, non-trivial solutions of a subcritical system in the presence of 
multiplicative noise can thus be destabilised through a "large" fluctuation. The 
trivial solution then becomes permanent. This simpler model explains the behaviour 
reported above in our dynamo simulations. 



Fig. 11. Potential G{X) as a function of the order parameter X. The potential represented here 
takes the form JSj setting = 5 and the prefactor of the term to 2 (solid line), 2.3435 (dotted), 
2.7 (dashed) et 3.7 (dot-dashed). 

The presence of a multiplicative noise implies that previously well defined equi- 
librium values (in the absence of noise) become probabilities. Let us consider a 
simple supercritical example of normal form 



If one assumes that ^(i) is a white and gaussian noise of variance Z?, one can then 
show that the equation governing the evolution of the probability P of X to have a 
given value is 



G(X) 




— = [^, + Cit)]X-X\ 



(6) 



dP 
'dt 



d 

dX 



(7) 



Integrating this equation in X for steady probabilities, yields 
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The right-hand-side vanishes for X = 0, which imphes that J = 0. The change of 
variable q — XP provides 

and thus P = C|X|2W^)-ie-^'/^. (9) 

uX D V X J 

It is worth noting that starting from a muhiphcative gaussian white noise (thus 
with a symmetric probabiUty distribution), one obtains a non-symmetric probability 
distribution for the variable X. 

One can show in a similar manner that for a subcritical bifurcation with a normal 
form given by ([4][5]), the probability P{X) to meet a given value is 

P^C^X^2(^./D)~1^xVD-X-/2D ^ (10) 

This simple example can, off course, not account for the full complexity of the 
coupled non-linear MHD system, yet it demonstrates how the structure of the in- 
duction equation accounts for the tendency for a stronger asymmetry in PDFs of 
the magnetic than in the kinetic energy (see figure 12). 




Kinetic Energy (^^ Magnetic Energy 

Fig. 12. Probability Density Functions for the kinetic energy (left) and magnetic energy (right), 
for Pm = 3 and E = S.IQ-*, Ra = 3.7 X Rac (top), E = 10"' *, Ra = 3.3 x Rac (bottom). The top 
line corresponds to a subcritical bifurcation, the bottom one to a supercritical bifurcation. 

2.5. Multiple Solutions 

In the course of the investigation reported here, we did not encounter proper multi- 
ple co-existing branches of dynamo action (but see section 3). We however observed 
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an interesting behaviour which happens to be directly related to the non-uniqueness 
of the purely hydrodynamic solution itself. Both the magnetic and the kinetic en- 
ergies are presented for E = 3.10^^, Pm = 6 and Ra = 2.57 x Rac in figure fT3l 
Two different simulations, which only differ by a minute initial fluctuation, are pre- 
sented on this figure. The behaviour is rather unexpected and untypical. While the 
simulations are first similar, the dynamo solution appears to reach a stable state 
(damped oscillations) but then violently destabilises. In one case dynamo action is 
maintained (left), but it is lost in the other one (right). 




Fig. 13. Magnetic (solid) and kinetic (dashed) energies versus time for E = 3.10^*, Pm = 6 and 
Ra = 2.57 X Rac. Left and right graphs correspond to two independent simulations with initial 
conditions which differ only from a very small perturbation (identified in the text as "run 1" and 
"run 2"). 



For this parameter space, we identified at least two stable solutions to the purely 
hydrodynamic problem. Further more, we observed (figure [T4|) . a modification in 
time of the dominant aziniuthal Fourier mode both for the kinetic and magnetic 
fields. Simulations were initialised with a random perturbation of the magnetic and 
temperature fields and a fluid at rest. The chosen parameters are here extremely 
close to the threshold value, and the magnetic field growth is slow (the growthrate 
is 0.058). Thermal convection therefore had enough time to develop largely unin- 
fluenced by the Lorentz force. For this parameter regime the kinetic energy of the 
most unstable mode of hydrodynamic convection is steady in time and the domi- 
nant modes are m = 6 and m — 0. When the field reaches saturation is has the 
same structure as the flow, dominated by an to = 6 symmetry. This mode however 
quickly destabilises and, as illustrated on figure I14| after a competition between 
the TO = 4, 5 and 6 modes, the dominant modes becomes (both for the kinetic and 
the magnetic energy) to = 5 for the first run and m = 4 for the second one. The 
simulation "run 2" , corresponds to a loss of dynamo action, as the to = 4 convection 
mode does not appear able to sustain the field against diffusion. Both simulations 
("runl" and "run 2") were carried further in time for another 30 t^, to assess the 
stability of their final state. 

This intricate behaviour therefore comes as the consequence of the existence of 
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three possible modes of hydrodynamic convection (all stables). The fastest growing 
mode, m — Q does act as a kinematic dynamo, but this dynamo is non-linearly 
unstable and destabilises either to the ?7i = 4 or to the m = 5-mode. The later is 
also a dynamo, but contrary to m = 6 it is a stable solution. The m = 4-mode 
does not act as dynamo, the magnetic field is then lost. We did not observe any 
heteroclynic cycle. The behavior reported here is reminiscent of self-killing dynamos 
as described in [TO] . 



3. Further developments 

We presented here a study of the onset of dynamo action over a wide parameter 
range. We should note however that one parameter was, for simplicity, kept constant 
throughout the entire study: the Prandtl number (Pr = v / k). It would be interesting 
to know how this parameter affects the picture presented here. 

Since this work was performed and published in french in 2005 [6] , we have fur- 
ther investigated the parameter space towards large values of the magnetic Prandtl 
number. Yet another kind of bifurcation diagram has been produced. It is character- 
ized by two co-existing branches of different energies. These results will be presented 
and discussed elsewhere [11] . 
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